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I review constraints on solar system-bound dark matter, and discuss the possibility that dark 
matter could be gravitationally bound to the earth and other planets. I briefly survey various 
■ empirical constraints on such planet-bound dark matter, and discuss effects it could produce 

p-^ ' if present, including anomalous planetary heating and flyby velocity changes. 

6 

+i , Cosmology suggests that only around 4% of the mass-energy density of the universe is ordinary 

I baryonic matter. In addition to ordinary matter, around 23% of the mass-energy density consists 

of gravitationally attractive "dark matter", and the remaining 73% consists of gravitationally 
repulsive "dark energy" . 

^ I Little is known about dark matter, other than that it is electrically neutral - hence dark, 

since it does not radiate photons - and that it interacts gravitationally through its mass-energy 
00 ■ density. Among the unanswered questions are: Is dark matter bosonic or fermionic? Is it self- 

I annihilating (i.e., its own antiparticle, or having equal abundances of particles and antiparticles) . 

' What is the mass (or are the masses, if more than one constituent) of dark matter particles? 

What are the non-gravitational interactions of dark matter? 
i Dark matter can be gravitationally bound on different scales. We enumerate possibilities in 

the following subsections. 



1 Gravitationally Bound Dark Matter 



1.1 Galactic Halo Dark Matter 

Our galaxy is believed to be surrounded by a dark matter halo, with a mass density (to within 
a factor of about two) of p ~ 0.3GeV/c^cm~'^. The galactic dark matter is believed to have a 
roughly Maxwellian velocity distribution of the form 

/a?;2g-3^,V(2^)2)^ ~ 270km s'^ . (1) 

In this formula, v = Wdm + "W© ; with Wdm the velocity of the galactic dark matter relative to the 
earth, and with the velocity of the earth through the galaxy. The latter has two components, 
the earth's orbital velocity (to which, at a low level, one should add the surface rotational 
velocity), and the velocity of the solar system barycenter through the galaxy. During part of 
the year these velocities add, and during part they subtract, leading to an annual modulation of 
the dark matter detection rate observed from earth which was proposed by Drukier, Freese, and 
Spergell^ and Freese, Frieman, and GouldEI as a signal for searching for dark matter. Possible 
observation of galactic dark matter through this signal has been reported by the DAMA/LIBRA 
experiment^. 



1.2 Solar System-Bound Dark Matter? 



In addition to dark matter bound to the galactic center of mass, it is possible that dark matter 
could be bound on smaller scales, such as to the sun or to planets. Dark matter bound to the sun 
would change planetary orbits, by giving additional perihelion precessions beyond that predicted 
by general relativity, and by changing the Kepler's third law relation between the period and 
semi-major axis of orbits. A number of authors (see Frere, Ling and Vertongen^^lj Sereno and 
Jetzer^, lorio^, and Khriplovich and Pitjeva^ have studied these effects, with the conclusion 
that the density of sun-bound dark matter is constrained by p < 10^GeV/c^cm~^. Because the 
earth's rotational velocity can add to the earth's orbital velocity, or subtract from it, depending 
on the sidereal time, sun-bound dark matter at densities well above the galactic density could 
be detected through a search for a daily sidereal time modulation in the DAMA/LIBRA and 
similar experiments. Such a modulation would have a 24 hour period, which would allow it to be 
distinguished from detector channelling effects of galactic halo dark matter, which as discussed 
in Avignone, Creswick and Nussinovi^ would lead to a sidereal time modulation with a 12 hour 
period. 



1.3 Earth-Bound (or Planet-Bound) Dark Matter? 

Another possibility is that there may be dark matter gravitationally bound to the earth. As 
shown in Adler^, one can place a direct bound on this, by using the fact that for a satellite of 
negligible mass in a circular orbit around a body of mass M , tracking to give the orbit radius 
R and orbit period T give the product GM of the Newton gravitational constant and the mass, 

GM = ^'k^B^IT'^ . (2) 

This formula, and its generalizations to include elliptical orbits and perturbations on an inverse 
square law force, can be applied as follows to give a useful bound on dark matter orbiting the 
earth. Consider first the LAGEOS geodetic satellite, a very dense nearly spherical satellite in 
orbit at a radius R ~ 12, 300km. Accurate tracking of LAGEOS, and modelling of drag effects, 
gives an accurate value for GM^, where is the mass of the earth plus the mass of any dark 
matter below the LAGEOS orbit. Similarly, accurate tracking of lunar orbiters gives GM^, with 
Mfn the mass of the moon. A more accurate way of getting the moon's mass (times G) is through 
tracking of the Eros asteroid flyby, with gives a very accurate figure for R(^/m = GM^/GMm- 
Finally, lunar laser ranging, which measures the relative distance between the earth and the 
moon, gives the product GMcombined for the earth-moon system, including the mass GAfdm of 
dark matter in orbit between the LAGEOS orbit and the moon's orbit, 

GMcombincd = GMe + GM„, + GM^^ . (3) 
Using this, and the earth to moon mass ratio obtained from the Eros flyby, we get 

GMdm ~ GAfcombincd - GM(s - GM(s/R(s/m , (4) 

where I have used an approximately equals sign because there could be a roughly 1% correction 
coming from dark matter in the vicinity of the moon or the Eros asteroid (see^ for details). 
Using the best current numbers (supplied to me by Slava Turyshev^^ one finds from equation 
(HI) that 

GMdm ~ (0.3 ± 4) X lO^^GMe , (5) 

with the dominant error coming from lunar ranging (which will be considerably improved over 
the next decade). 7f this bound were attained, and if the dark matter mass were uniformly 
distributed below the moon's orbit, the density would be /9 ~ 6 x 10^°GeV/c^cm~^, which is a 



much higher density than both the galactic halo dark matter density and the upper bound on 
the density of sun-bound dark matter. So in principle, there could be substantial amounts of 
dark matter gravitationally bound to the earth (or generalizing from this, to other planets). 

2 Possible Applications of Earth- and Planet-Bound Dark Matter 

What I now want to say is necessarily speculative, since it is conditioned on the possibility that 
there may be significant amounts of dark matter gravitationally bound to the earth and other 
planets. 

2.1 Jovian Planet Anomalies 
In AdlerE] 

I suggest that dark matter accretion may play a role in explaining certain features 
of the Jovian planets Jupiter, Saturn, Uranus, and Neptune. From de Pater and Lissauer one 
learns that the surface heat fluxes H of these planets, in ergs/cm^s, are 5440 (Jupiter), 2010 
(Saturn), <42 (Uranus), and 433 (Neptune). The larger heat fluxes have not been completely 
accounted for (for example, through internal radioactive decays) and it is possible that a new 
mechanism may be at work. I suggest that the accretion of planet-bound dark matter may 

account for unexplained internal heat production. This would be possible if each planet has an 
associated dark matter cloud (perhaps linked, as suggested in^, to formation of the planet), and 
if the dark matter is accreted with a low energy release efficiency, so that bounds on earth heat 
production arising from accretion of galactic dark matter are not violated. Another peculiar 
feature of the Jovian planets is that Uranus, which is otherwise very similar to Neptune, has 
a much lower internal heat production, as well as being the only planet which has its rotation 
axis lying on its side rather than nearly normal to the ecliptic. I suggest that the collision that 
is hypothesized to have tilted the axis of Uranus could also have knocked Uranus out of its 
associated dark matter cloud, accounting also for the anomalously low internal heat production. 

2.2 Flyhy Anomaly 

When spacecraft are put into near-earth hyperbolic orbit flyby trajectories, there is a several 
hour segment when the spacecraft is close to the earth and cannot be tracked. When the 
incoming orbit is extrapolated, using the best orbit fitting programs, to predict the outgoing 
orbit, a discrepancy between the predicted and observed outgoing velocity is found, as reported 
in Anderson et al.^^^For example, the Galileo II flyby on 12/8/92 shows a velocity discrepancy 
of -4.6 mm/s (with an estimated error of 1 mm/s), that is the outgoing velocity is lower than 
expected, while the Near flyby on 1/23/98 shows a velocity discrepancy of 13.46 mm/s (with 
an estimated error of only .01 mm/s), that is the outgoing velocity is higher than expected. 
The discrepancies are roughly of order 10~® of the total velocity, so constitute a relatively large 
effect, and as we have just seen, can have either sign. 

There are at least four possibilities for explaining these flyby anomalies: 

i. The effect is an artifact, resulting from the omission from the orbit fitting programs of known 
physics. However, no one has been able to identify a candidate so far, since all the obvious things, 
including moon and planetary perturbations, special and general relativistic corrections, solar 
wind effects, Van Allen belt collisions, tides, and thermal effects, have been taken into account, 
or estimated and shown to be too small to be relevant. 

ii. New electromagnetic physics is involved. However, this seems unlikely, since experimental 
tests of quantum electrodynamics and vacuum linearity have a much higher accuracy than 1 part 
in 10^. 

iii. New gravitational physics is involved. This would have to be quite unconventional. MOND 
(modified Newtonian dynamics) , which suggests a modification in Newtonian gravitation for very 



low accelerations, does not predict the flyby anomalies. I (and others) have taken a look at the 
parameterized post-Newtonian (PPN) formalism that gives a phenomenology for studying (in 
the weak gravitational field, low velocity regime) metric theories of gravitation that differ from 
general relativity. Using the modified equation of motion given in WiliE! and the 2006 bounds 
on the PPN parameters, I find that all terms are too small, by several orders of magnitude, 
to account for the flyby anomaly. This still leaves t he p ossibilitv of non-metric gravitational 
theories, and some ideas for these have been proposed^"'^-l 

iv. The effect comes from collisions of the fiybys with earth-bound dark matter This 
possibility is analyzed in my paper AdlertEl. Velocity decreases can arise from ordinary drag, 
as a result of elastic scattering D{mD) + N ^ D(mD) + N, with m£) the mass of a dark matter 
particle D and a spacecraft nucleon. Velocity increases can arise from exothermic inelastic 
scattering D{m£,) + N ^ D'{m£)i) + N, with the primary dark matter particle of mass rriD 
scattering into a secondary of mass < m£,, and with the secondary recoiling predominantly 
backward, so that the nucleon gets a forward kick. 

My analysis of the dark matter possibility shows that the following strong constraints are 
required, in order for dark matter collisions to provide an explanation for the flyby anomaly: 

i. The dark matter in orbit around the earth must be localized well within the moon's orbit and 
not too near the earth. 

ii. The dark matter mass must be much smaller than a GeV. 

iii. The dark matter scattering cross section on nucleons must be high, of order 10^'^'^ to 

iv. The dark matter must be non-self-annihilating and stable in the absence of nucleons. 
These constraints require a form of dark matter that is very different from the customary 

assumption that dark matter is a multi-GeV lightest supersymmetric particle. 

My current work consists of modelling dark matter orbiting the earth, with two species 
in different configurations, one of which elastically scatters from nucleons, and one of which 
inelastically scatters, to try to fit the Anderson group data. I assume dark matter distributions 
obtained by averaging a circular orbit of radius r at an angle x relative to the earth's rotation 
axis, over precession around the earth's axis, giving a truncated spherical shell of radius r, height 
2|z| = 2r sinx, and relative density l/-y/(r sinx)"^ — z'^- For a given point characterized by r, z, x 
there are two relevant dark matter velocities that correspond respectively to an up-moving and 
a down-moving dark matter orbit segment passing through the point, and the corresponding 
velocity vectors are easily worked out. From these velocities and the spacecraft velocity, one 
can calculate the cross section-averaged velocity change imparted to a spacecraft nucleon, and 
integrating the corresponding work increment over the spacecraft orbit gives the final energy 
and velocity change. I currently have the core programs for this calculation written, and now 
have to do a systematic search of the parameter space to try to find a reasonable fit to the data. 
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